A Non-equilibrium STM model for Kondo Resonance on surface 
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Based on a no-equilibrium STM model, we study Kondo resonance on a surface by self-consistent 
calculations. The shapes of tunneling spectra are dependent on the energy range of tunneling 
electrons. Our results show that both energy-cutoff and energy-window of tunneling electrons have 
significant influence on the shapes of tunneling spectra. If no energy-cutoff is used, the Kondo 
resonances in tunneling spectrum are peaks with the same shapes in the density of state of absorbed 
magnetic atoms. This is just the prediction of Tersoff theory. If we use an energy cutoff to remove 
high-energy electrons, a dip structure will modulate the Kondo resonance peak in the tunneling 
spectrum. The real shape of Kondo peak is the mixing of the peak and dip, the so-called Fano line 
shape. The method of self-consistent non-equilibrium matrix Green function is discussed in details. 

PACS numbers: 68.37.Ef, 73.23.-b, 72.10.-d, 73.20.-r 



I. INTRODUCTION 

Scanning Tunneling Microscope (STM) has been 
widely used to study atomic and electronic structures 
of surface and atoms absorbed on it. A currently inter- 
esting topic is Kondo effect of magnetic atoms on non- 
magnetic metal surface. The interest along this direction 
is motivated by recent advance in the construction of 
novel nano-structure device on surface. STM has been 
used to probe the Kondo resonance of a single magnetic 
atom.pl H H II IE H The magnetic cobalt phthalo- 
cyanine molecule absorbed on Au(lll) surface shows 
Kondo effect by having cut away eight hydrogen from the 
molecule. 7] Kondo resonance haven been found in other 
magnetic molecules absorbed on metal surfaces such as 
Co(CO)„ molecules on Cu(100) surface ||. Additionally, 
STM has also been used to study the formation of Kondo 
molecule. 0, ^| The most interesting phenomenon is so 
called Quantum Mirage 0, 0] due to the refocus of 
Kondo resonance on surface. 

Generally STM tunneling spectrum isn't exactly cor- 
respondent to the atomic and electronic structure of the 
surface. To obtain the realistic atomic and electronic 
structure from STM tunneling spectrum, it is valuable 
to establish a believable STM theoretical model that is 
related the realistic atomic structure or electronic struc- 
ture to the STM tunneling spectrum. Based on the Ter- 
soff alike theory, [TJl [Tj, [TjJ [TJJ the differential conduc- 
tance of electronic tunneling in STM is proportional to 
the Local Density of State(LDOS) of surface state. The 
resonance peak at zero-bias voltage in dl/dV curve is 
generally the mixing of Kondo resonance with other reso- 
nance such as d-resonance from 3d magnetic atom. Q The 
STM models based on Anderson Model have successfully 
explained the Fano shape of Kondo resonance using the 
equilibrium many-body theory 0, 0, 0, [2jJ, EI and 
non-equilibrium theory. [22| 
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Quantum tunneling of electrons in STM device is out 
of equilibrium in nature. In this paper, using a non- 
equilibrium self-consistent method, we calculate the STM 
tunneling spectra for magnetic atoms on metal surface. 
Kondo effect can be solved using analytical or numer- 
ical methods, such as, the Non- Crossing Approxima- 
tions (NCA) and other Larger N methods are suited for 
the infinite U Anderson model. [23J, |2J, |2£j The numeri- 
cal renormalization group (NRG) and the exact Bethe 
Ansatz need to linearize the dispersion relation near 
Fermi energy after having simplified into one-dimension 
model. [2l| The method used in this paper is similar to the 
iterative perturbation theory (IPT).[2E E3 The big dif- 
ference is that in this work the Kondo resonance peak is 
decoupled from other parts of the Green function. How- 
ever for IPT, the Kondo resonance is imported into the 
Green function of magnetic atom by interpolating Sq self- 
energy into an appropriate self-energy ansatz. The de- 
coupled scheme in this work can be obtained by more 
detail theoretical analysis such as in reference. |2q 

As a company of Kondo resonance, the so-called zero- 
bias anomaly has contribution to the differential con- 
ductance near zero-bias voltage. The zero-bias volt- 
age anomaly may have different origins such as the or- 
thogonality catastrophe 29] and the co-tunneling mech- 
anism [3(1 HSI ■ The electrons have no enough energies 
for most of many-body events at small bias voltage. The 
zero-bias anomaly appears as a local minimum in dl/dV 
curves near zero-bias voltage. 

Experimentally Kondo resonance in dl/dV curves 
are clear dip structures for single absorbed magnetic 
atoms pj 0, 0, butpeak structures for absorbed mag- 
netic molecules 0, [g- In this paper, the signals of 
Kondo resonance from the dl/dV curves are Fano-line 
shapes and are the mixing of Kondo resonance and dip 
structure due to zero-bias anomaly. This is well con- 
sistent with experimental observations. The zero-bias 
anomaly generally appears as a dip structure near zero- 
bias and is frequently covered by background noise. If 
other resonance peak appears at zero bias, the zero-bias 
anomaly can shows itself by modulating the resonance 



2 




FIG. 1: The energy-level distribution of the STM device. 
Two horizontal lines labelled with fi T and n s illustrate the 
chemical potentials of the tip and the surface. The fi s is the 
chemical potential for the surface where the Kondo resonance 
happens. The red areas show the states of the tip and sur- 
face, which are filled with electrons. The chemical potential 
fi s of the surface is fixed to zero. The difference of the chem- 
ical potentials for the tip and the surface is induced by the 
bias voltage V& applied to the surface. Two thick bold lines 
present the splitting energy-levels for magnetic atom due to 
the on-site coulomb interaction. Two curved dash lines show 
a cotunneling event in the STM device. 



peak. This is why Kondo resonance generally appears 
company with a local minimum in dl/dV curves. If we 
properly reduce the back ground noise, the dip struc- 
ture also can be observed in the tunneling spectra. The 
shapes of tunneling spectra are dependent on the energy- 
range of tuneling electrons. The paper is organized as 
following, after brief introduction we present the theo- 
retical STM model in section two. In the section three, 
we present the non-equilibrium self-consistent method in 
detail. Our main results are presented in section four. 
Two appendixes are at the end of the paper which in- 
clude the brief introduction of the matrix representation 
of non-equilibrium Green function and the derivation of 
the tunneling-current formulations. 



II. THEORETICAL STM MODEL 

The Non-equilibrium Green-function method 
(Schwinger-Keldysh Contour) [H |H |M 13 
has been used successfully to study the mesoscopic 
quantum transport of quantum dots and others meso- 
scopic quantum device. |37l 1381 l39l Phi Hll |42| The 
model Hamiltonian is splitting two parts H = Hq + Ht- 
Ht is the tunneling Hamiltonian which expresses the 
electronic transport between different parts of the STM 
device such as tip, surface and absorbed atom shown in 

[Fig.n. 
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p',P, d^ and p, t, d are the electron creation and anni- 
hilation operators of surface, tip and absorbed atom re- 
spectively. Vnaa-', B nk(J(J > and A kaa , are the tunneling 
matrix elements or hybridizing matrix elements between 
the tip and absorbed atom, the tip and surface, and the 
absorbed atom and surface, k is the energy-level index 
of the surface, n the STM tip and a is the spin index 
of electron. e pk(J , £tna and e a are the energy levels of 
the surface, the tip and the absorbed atom respectively. 
U is the on-site Coulomb energy of the absorbed atom. 
The hybridization Hamiltonian, the second term of Ht, 
has the same form as the tunneling Hamiltonian. Thus 
hybridization can treat the same foot as the tunneling 
events. 

The STM probes the tunneling current or the change 
of charges in the tip when a bias voltage applies to the 
surface. The operator of electron number at the tip is 
written as Nt — VL„ t and the electronic current is 

' l — /na na n<T 

its time-derivative of J T = -e^f = ^[N T (t), H], After 
having been calculated the commutations the current is 
expressed as 

{J T ) = %Re(EV nnl (tld w ,)<+ ]T B n ^ r (tlp ka ,)<U3) 

naa' nkaa 1 

where the lesser Green function G < (t, t') — i < 
A(t'Y[B(t) >. The similar formula has appeared in 
Ref.[22j but replaced the real part with imaginary part. 
Two formulas are the same if having considered without 
image unit i in the definition of lesser Green function in 
that paper. In the appendix B, we derive bellow current 
formula using the matrix non-equilibrium Green func- 
tions. 

(J T > = J1 + J2 (4) 
J x = ^ f dwi?e7>[G ct(r) H£ ct(<) H] 

J c 

J 2 = ^^dwi?eTr[G ct(<) (w)S ct(Q) (w)]. 

where Tr presents the summation for all energy-level and 
spin indexes J2nn'aa'- ^ ne superscript c indicates that 
Green functions are matrixes defined in the appendix A. 
The r, f , < and > indicate the four elements of the 
matrix Green function Eq. (|A1|) and the r and a are the 
retard and advance Green function defined in Eq. I|A4|) . 
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equations for the tip, adatom and surface, which are ex- 
pressed as 



Q, lbn G ct = ff ct M + ff ct Mi] ct MG ct M (5) 



V V* B B * 



7 cd = g cd (w)+g c ' i (w)S c ' i (w)G cd ( 
G cp = g cp (uj) + g cp (uj)Z cp {oj)G cp (uj) 



V* A B* 



or diagrammatically illustrated in FigEl We replace the 
C-i (c) 2 free Green functions in the self-energy formulas Eq. I|B5|) 



(a) ° %"""?, mmmMmmm with new Green functions 



FIG. 2: (a) Feynman diagram representation of Dyson equa- 
tion for the tip. The Dyson equations for absorbed magnetic 
atom and surface are similar, (b) The physical event for the 
first term of the tip self-energy £2 includes two tunneling 
events. The electron jumps to the absorbed atom and re- 
turns to the tip. (c) The physical events for the first term of 
the tip-energy E3 includes a cotunneling event from the tip 
to the surface via the absorbed atom and a single tunneling 
to return to the tip from the surface. 



The current Ji flows into the tip from the surface and 
absorbed atom, and J2 flow out of the tip to the surface 
and absorbed atom. The limitations of integral of the 
current is dependent on the problem we will consider and 
can choose Eq. (|B7|) or Eq. (|B8() . Near zero bias we need 
an energy cutoff and Eq. (|B8() is a good choice. 



III. THE METHODS OF NUMERICAL 
CALCULATION 

The non-equilibrium Green function is defined on a 
contour with two branches: the positive branch from 
negative infinite on time axis to positive infinite and the 
negative branch from positive infinite return to negative 
infinite. The contour green functions can be replaced by 
a matrix Greens functions. The matrix Green can be 
calculated using the standard perturbed methods in the 
same manner as the zero-temperature Green function in 
the standard textbook. j3^ In Appendix A, we present 
the main formulas and their computational methods. 



A. The key formulas for self-consistent calculation 

We assume the Hamilton Ho can be solved exactly and 
the tunneling Hamilton treats as perturbations. The self- 
consistent calculations are used to obtain the full Green 
functions. At first, we solve the Hamilton Ho and the ob- 
tained Green functions are used as the free Green func- 
tions for self-consistent calculations. At the next step, 
the self energies are calculated using the Eq. I|B5|I . Three 
new Green Functions are calculated using three Dyson 



n' n" a' a" 

+ Y A. , , A*„ „ G cp H + EgUw) 

k' k" a' a" 

£*7 aM = V B , , B* G ct (u) 

n'n"{T' a" 
a'a" 
a'fj' 
k<l<a<f3< 

The Eq. 10 and three Dyson equations Eq. JSJ) form a 
self-consistent loop. In Fig[3and Eq. (|B6|> . the S3 self- 
energies represent the cotunneling mechanism which will 
contribute part of the zero-bias anomaly in dl/dV curve. 

B. The free Green functions for the surface and tip 

At first, we calculate the equilibrium retarded and ad- 
vanced Green functions as the free Green functions used 
in following non-equilibrium self-consistent calculations. 
The surface is treated as two-dimensional electronic gas. 
The density of state (DOS) of two-dimensional electronic 
gas is flat and featherless with 4.0 eV width centered at 
Fermion energy of the surface. The well-shaped sharp 
STM tip is constructed with smaller number of atoms. 
The tip can be modeled as an atomic cluster with dis- 
crete energy-levels. Because the electron number is lim- 
ited, there exists a highest-occupied energy-level such as 
the HOMO (Highest Occupied Molecular Orbit) energy 
E HOMO for molecular systems or Fermi-energy for metal, 
all states with energies less than the highest-occupied 
energy-level are occupied, and all states with energies 
larger than the energy-level are empty. The electrons 
with energies near the energy-level are important to the 
physical and chemical properties of atomic clusters. Be- 
cause the tip is in the metal environment, we assume the 
tip is in metal state and has a sharp Fermi surface. The 
tip models as a single energy level identified as Fermi en- 
ergy [It- The free retarded Green function for the tip is 



4 



written as g^ a = UJ _ £t ^ +iSt , where e tna = \i T . The DOS 
of the tip has the Lorentz shape and centers at its Fermi 
energy pt with width S t . 

C. The free Green functions for magnetic atom 



TABLE I: The parameters for numerical calculations (unit 
eV). The band width for surface state is 4.0 eV 



A 


V 


B 


u 


£d 


5 d 


St 


0.20 


0.10 


0.05 


2.80 


-0.80 


0.30 


1.00 



The main term produced the Kondo effect is the spin- 
flipping term, that is the last term in the Hamilton 
Hq. We assume the Hamilton had been solved exactly. 
The free Green function of magnetic atom including the 
Kondo resonance is calculated using the Hartree-Fock ap- 
proximation with U 2 and U 4 order self energy. ^3 The 
retarded Green function of the Anderson impurity is writ- 



ten as G r da (tjj) = LU _ E d +iA where e d = e + n s U. From 
the results of particle-hole symmetric Anderson model in 
reference, 0, El up to the U 2 and U 4 self energy, the 
Gk (w) near zero frequency is written as 



and the surface, and treated as an input parameter. The 
scheme Eq. J7J) isn't simple decoupling because the width 
and height of the decoupling Kondo resonance are de- 
pendent on the energy-level width and splitting of mag- 
netic atom. The Friedel sum-rule Eq. also tights the 
Kondo resonance and other parts of Green function to- 
gether in the following self-consistent calculations. At 
present, the scheme only used in the study of Anderson 
model in Kondo regime. 



D. Summary of self-consistent calculations 



(7) 



E 2 (o;) - -a^-z a/ A((^) 2 + (^) 2 ). 



where a R = (3 - ^)(^) 2 + 0.0553(^) 4 , a 7 = 
\ [(^) 2 + 6(5- f)(^) 4 ] and A = np Q \A\ 2 . We find 
from our calculations that it is enough up to U 2 term. 
At zero temperature, the height of the peak is ex- 



actly po 



£2(0;) satisfies the general Fermi-Liquid 



rfJmS 2 (0) 
du 



at zero 



properties im£ 2 (0) = and 

temperature. It is also a good approximation to asym- 
metric Anderson model in Kondo regime if the asym- 
metric Anderson model is still described approximately 
by Fermi liquid. The approximation has been used to 
study Kondo resonance on surface. The main temper- 
ature effect is introduced by Eq. (J7J and other tempera- 
ture effects are introduced by Fermi distribution function 
n f (uj). 

The Green function of magnetic atom, similar to the 
reference, 20 is written as 



gl,(w)=§ r da (w) + Z t ,G K (w-6E), 



(8) 



where Gk{u) is Kondo resonance peak and SE is the 
shift of Kondo peak away from the Fermi energy. SE 
is generally very small about several meV. The choice 
of Z a , (Z-[,Zi)={1.0, 0.0), meets the condition that the 
height of Kondo peak is equal to ^ in the total den sity 
of state of magnetic atom. Instead of in reference |2Cj 
where the Kondo resonance is Lorentz shape, in this 
work the treatment of Kondo resonance is improved with 

non-Lorentz shape. The Green function q 1 ! = } , ... , 

where Sd doesn't include the contribution from the hy- 
bridizations with surface state and tip and doesn't con- 
tribute to Kondo resonance. 8 d can be seen as the con- 
tributions from the bulk electrons j^TJ, [5(j or other com- 
plex interactions between the absorbed magnetic atoms 



Up to now, we have known the free retard Green func- 
tions for the tip, surface and magnetic atom. Based 
on Eq. QA5J1 . we can transform the retard Greens to 
the Matrix Green functions as the free matrix Green 
functions for non-equilibrium self-consistent loops. In 
this work, the trial Green functions Gq(li) for the self- 
consistent calculations are equal to their free Green func- 
tions Gq(u>) — g c (u). In the self-consistent loop we cal- 
culate the occupation number 



+00 



duf(w - (is)ImG da (uj). 



(9) 



the new position of energy-level of magnetic atom e d a — 
e d + n s U and the energy-level width A = —ImY. d (uo = 
Us) which is the sum contributed from the hybridizations 
with surface state and tip. Thus the shape of Kondo res- 
onance is also changed in each loop. If the changes of 
occupation numbers are smaller than 0.001, the occupa- 
tion numbers, Green functions and self-energies converge 
simultaneously. The Friedel sum-rule Eq. is satisfied 
and the self-consistent calculations stop. 

The spectral functions are obtained by formula A(u>) = 
\{G < {lo) — G > {lo)) and the density of State p(u>) — 
^A{ui). The tunneling current is calculated using Eq. 
The differential conductance is numerically calculated 
from the curve of the tunneling current via bias volt- 
age. In our self-consistent perturbed calculations the ba- 
sic equation G T + G T = G > + G < for the matrix Green 
function is only approximately satisfied with accuracy 
about from 0.001/eV to 0.01/eV. 



IV. 



RESULTS AND DISCUSSIONS 



In this work we concentrate on the small bias volt- 
age. It is convenient to drop the indexes of V naa i , B n k aa i 
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FIG. 3: The density of state for magnetic atom (a,c) and 
the surface (b) with A=0.20 eV and U=2.8 eV. The broad 
peak in (a) is the wide energy level for magnetic atom; the 
small sharp peak is the Kondo resonance. The solid line is the 
density of state (DOS) at 40 K and dot line at 240 K. The 
surface density of states and the detail structure of Kondo 
resonance are shown in the insert figure. 
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FIG. 4: (a). The tunneling spectra of magnetic atom con- 
stant energy-window [-2.0 eV , 2.0 eV] at near zero-bias volt- 
age at temperatures 40K and 240 K. Kondo peak becomes flat 
at 240K above Kondo temperature about 100 K. The tunnel- 
ing spectrum of empty surface at 40 K also plots in this figure, 
(b) The tunneling current at 40 K and 240 K when the tip on 
top of the magnetic atom. 



Calculations of Density of state (DOS) and 
Differential conductance 



The density of state of the absorbed magnetic atom 
and surface are shown in Fig. [3J From this figure [Fig. [31 
(a)], we can find one broad peak of energy level of 
the magnetic atom. The sharp peak at Fermi energy 
of the surface is well known Kondo resonance. The 
Fig. I^Jc) shows the flat effects of Kondo resonance, 
when temperature increasing above the Kondo temper- 
ature Tk- We estimate the Kondo temperature based 
on the width of Kondo Peak at zero temperature and 
find T^-=100 K. The Kondo temperature Tk also is 
estimated based on the equation Dyf7pe~~ 0, |44| . 
where p is the density of state of surface at Fermi energy 
J = A/(np)[l/\e d \ + l/\e d + U\] obtained by Schrieffer- 
Wolff transformation [45| and D is the effective band- 
width of surface state contributed to Kondo effects. We 
estimate the Kondo temperature T#-=123 K if 2D = 4 
eV. The 2D is approximate to the energy windows in our 
calculations from -2 eV to 2 eV. The Fig. [31(b) shows the 
density of state of the surface which is almost flat only 
there is a valley at the Fermi energy. The depression of 
surface DOS near Fermi energy is because some surface 
electrons at Fermi energy bind with the Anderson im- 
purities and screen its local spin by forming spin-singlct 
binding state. 

The peak structures of Kondo resonance have been 
found in the tunneling spectra of cobalt atom embedded 
in larger molecules absorbed on metal surfaces. 0,0- We 
calculate the differential conductance using the parame- 
ters in table JIJ) at temperature 40 K. The energy-range 
in the integral is from -2.0 eV to 2.0 eV, approximately 
from — oo to +oo in our calculations. We also calculate 
the tunneling spectrum of empty surface. From Fig. 0] 
(a) we can see that a peak appears at zero bias and it 
became flat at 240 K. It's just Knodo resonance in tun- 
neling spectrum. If the tip move away from the magnetic 
atom to empty parts of the surface, the tunneling spec- 
trum is flat. The peak structure of Kondo resonance 
and the flat spectrum for empty surface are celebrated 
because dl/dV curve is correspondent to the density of 
state and close to the prediction of the standard Tersoff 's 
STM theory [HIT^. 



Influence of energy-window and energy cutoff 
on shape of resonance peak 



and Akaa' and simplify to V, B and A without introduc- 
ing confusion. The model parameters arc collected in 
table JIJ). The parameters are close to Kondo regime 
(U >■ A < Ed) and close to the values for Co atom 
on Au(lll) surface provided in reference. pfj| . In Kondo 
regime, there is sharp Kondo resonance peak at Fermi 
energy besides the bare virtual resonance peak in the 
density of state of the absorbed magnetic atom. 



The dip structures are more common in STM tunnel- 
ing spectrum. The energy-window of tunneling electrons 
has important influence on the shapes of dl/dV curves. 
Several theoretical works [lg, |2l| have shown that, in 
order to obtain correct tunneling spectra, the calcula- 
tions of tunneling matrix need to use an energy-cutoff 
above Fermi energy to remove high-energy surface elec- 
trons. This is implicit that only part of electrons with 
proper energy-range can tunnel into the tip. At zero 
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FIG. 5: (a) The energy-window influences the shape of dl/dV 
curve near zero-bias. Three different energy- windows are used 
{-(eV b + SE),(eV b + 5E)] 5E=0.0, 0.0034, 0.02 (eV). (b)The 
dl/dV curves for single electron transmits through a square 
barrier with height Eb- The sharp dip structures are shown 
for small barriers. The insert figure illustrates the simplified 
STM device without absorbed atoms on surface. En = „ ■> 
is a characteristic energy for the barrier 



temperature, there exists clear energy-window with the 
Fermi energies of the tip and surface as its up-boundary 
and down-boundary. Only electrons with energies in the 
window can transmit to the other side of the barrier be- 
tween the tip and surface and tunnel into the tip. The tip 
can't capture high-energy electrons although they have 
already transmitted the barrier. 

At finite temperature, the energy distribution of elec- 
tron becomes smoothly across the Fermi energy. We may 
add a skip-energy SW to define energy-window at finite 
temperature to include more electrons and cutoff elec- 
trons with higher energy. We alter the energy-window 
to [-(eV b + SW),(eV b + SW)} and calculate the dl/dV 
curve again. The width of energy- window changes with 
the bias voltage. In fact cutoff effect of Fermi energy has 
already included in our current formula by lesser Green 
functions and self energies for non-interaction case. After 
self-consistent calculations the electronic interactions of 
different parts are introduced. The cutoff effect of lesser 
Green functions becomes weak even at very low tem- 
perature and Fermi distribution still allows high-energy 
electrons above Fermi energy. We use the energy cut- 
off eVb + SW to neglect electrons with higher energies. 
The integral limits are symmetric around Fermi energy. 
From Fig. 0(a), we can see that the dip structures appear 
when SW < 0.02 (eV). For SW = k B T = 0.00334(el/) < 
0.02(eF), there is a sharp dip structure. We can also 
find from Fig. EJ(a) that energy- window effect is only ef- 
fective near zero bias. If the minimum value of effective 
energy- window is larger than 0.02 (eV), the influence of 
energy- window can be ignored. 

Both the energy-window and energy-cutoff effects to- 



gether determine the local minimum of dl/dV curve at 
zero bias. We must analyze them more deeply. At large 
bias voltage V& the energies of most tunneling electrons 
are in energy-range from eV to eVf,. The number of 
electrons with high-energies is relative small. At very 
small bias voltage, the number of electrons within the 
energy-range is deceased to very small value. At this 
time, if we remove high-energy electrons the total dif- 
ference conductance will has a large decrease. Only SW 
less than 0.02 eV, in another word, more high-energy 
electrons are removed the local minimum near zero-bias 
is more prominent. 

We study a simplified STM model, that is, electrons 
transmits though a single barrier which is vacuum region 
between the tip and surface. From the insert figure of 
Fig- El (b), the height of the barrier is Eb and the energy 
of incident electron is E. By standard quantum mechanic 
calculation, the transmitting coefficient is 



T(E) 



1 



B 



sinh( 



2ma 2 
h 2 



(E B -E)) + l 



(10) 



for E < Eb where a is the width of the square bar- 
rier. The differential conductance is calculated by the 
Landauer-Buttiker type formula G = =f-T(eV&). From 
the Fig. |3] (b) . We find the sharp dip structure for small 
barrier. The transmission coefficients increase with ener- 
gies of tunneling electrons. Compared with high-energy 
electrons, near zero bias the energies of tunneling elec- 
trons are small and the transmission coefficients are also 
small. 

The dl/dV curves obtained from our model have sim- 
ilar shapes to the dl/dV curves for the single barrier 
model is indicative that the dip structures are because 
the energies of tunneling near low bias locate in classic 
forbidden regime of the vacuum barrier between the tip 
and surface. We remove the high-energy electrons by in- 
troducing energy-cutoff is equivalent that we restrict the 
energies of tunneling electrons below the effective barrier 
and within classic forbidden regime near zero bias. The 
number of high-energy electrons is small, but their trans- 
mission coefficient is large. They will have important in- 
fluence on differential conductance because if at low bias 
most of electrons have low energies and small transmis- 
sion coefficient. Although STM device has simple struc- 
ture, many things are unclear such as the shape and size 
of tip. Generally the tip can't capture all tunneling elec- 
trons. At low bias, most of electrons locate in low-energy 
regime. The high-energy electrons and low-energy elec- 
trons are in different energy-scale so high-energy events 
are irrelevant to low-energy events. This is the reason 
why we introduce energy-cutoff to remove high-energy 
electron. The realistic STM devices are more complex. 
We must known the exact information of surface state. 
Of course, the directions of velocity of tunneling electrons 
also determine whether the tip can capture these tunnel- 
ing electrons. Thus we suggest that, to more close to 
realistic STM, we should add a captured probability fac- 
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FIG. 6: The differential Conductance vs bias Voltage at 40 
K with A=0.20 eV, U=2.8 eV. The solid lines with circles 
illustrate the results at 40 K. The Kondo temperatures Tk are 
about 100 K. The solid lines indicate the Fano shapes fitting 
the dl/dV curves at 40 K. The figure (b) is the asymmetric 
resonance when the Kondo peak shifts away the Fermi energy 
about 0.005 eV. (c) the dl/dV curve for nonmagnetic impurity 
and (d) for empty surface. 



tor W(ui) into the current formula, and it can be written 
as J = 3~ J W(uj)T(uj)duj, where T{ui) is transmission co- 
efficient. In this paper, W(uj) = 9(tu — n)0(fi + eVb — w). 
This is just equivalent to the introduction of energy- 
window in our calculations. This is also equivalent to 
limit all tunneling matrixes with non-zero values only 
in this energy- window just like in reference. [Tftl l2ll | Be- 
sides energy, W{...) could include position, momentum 
and any other more complex properties of STM device. 

Based on our results, the dip structures in dl/dV curve 
near zero-bias are induced by the tunneling among the 
tip, absorbed atom and surface because the energies of 
tunneling electrons are smaller than heights of effective 
barriers of corresponding tunneling process. Not only 
most of many-body events but also subsequent single tun- 
neling events are suppressed. The electronic tunneling 
events include subsequent single tunneling, co-tunneling 
and more complex tunneling process. If the background 
noises are reduced by proper method, the dip structure 
can be found in the dl/dV curve. If there are magnetic 
atoms on surface, Kondo resonance will enhance the dif- 
ferential conductance near zero bias, the dip structure 
has chance to be observed by modulating the shape of 
Kondo resonance. Thus the most of Kondo resonance 
observed in STM tunneling spectra are modulated by a 
local minimum. 



C. Calculations of differential conductance with 
proper energy-cutoff 

In this section we will study STM tunneling spectrum 
for the same parameters as in section A but we use a 
energy-cutoff from eV to eV 6 (see Fig.JTJ). The Eq. l|B8|) 
is well common used formula for the calculations of elec- 
tronic transport near zero temperature. The results are 
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—v— (2) No Kondo peak 

—a— (3) No L 3 and Kondo Peak 
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FIG. 7: (1) Open squares show the dl/dV curve for the full 
calculation including Kondo peak and (2) solid line with open 
down-triangles show the dl/dV curve without Kondo reso- 
nance ( Eq. 10) in self-consistent calculations. The (3) up- 
triangles for without both £3 self-energies and Kondo peak, 
and the (4) diamonds for empty surface are also plotted in 
this figure. The contribution of Kondo resonance to dl/dV 
curve (bold thick-solid line) is obtained by having subtracted 
(2) from (1). 



consistent with most of STM tunneling spectrum of mag- 
netic atoms on metal surface. The differential conduc- 
tance (dl/dV) is plotted in Fig|Hfa). After having fit- 
ted the dl/dV curves at 40 K with Fano shape ^2±£L 



1+e 2 



where s 



k B T q 



we get Fano parameter q=0.0 for (a) 



which is good fitting with Fano line-shapes. The symmet- 
ric resonant peak has been found in the STM tunneling 
spectra of Co/Au(lll) and Co/Ag(lll) Q. Experimen- 
tally, some STM tunneling spectra of magnetic atoms on 
fcc(lll) surface of novel metals are asymmetric dips 0,0 
or peaks. 0. We shift the position of Kondo resonance 
about SE = 0.005 eV away from the Fermi energy. We 
obtain asymmetric resonance peak very close the shape 
observed in experiments [lj, |3J and q values by fitted with 
the Fano-line-shape is 0.30 for FigHJb). 

In the nonmagnetic calculation the parameters are the 
same but {7 = 0. The clear dip structure is found from 
the dl/dV curve (FigOJc)). For the empty surface with- 
out absorbed magnetic atoms the dl/dV curve has also 
a dip structure near zero-bias. There is no dip structure 
in the density of the state of the nonmagnetic impurity 
and empty surface at the Fermi energies. Thus dip in 
dl/dV curve is the tunneling phenomenon such as single 
tunneling and co-tunneling events. The magnetic impu- 
rity and even nonmagnetic can enhance the dip structure. 
The resonant structures near zero bias for nonmagnetic 
impurities are found in other non-equilibrium calculation 
such as in reference. [22| 
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D. The detail analysis of dl/dV curves 



We analyze the dl/dV curves in detail. In the current 
formulation Eq. (@J, the tunneling current is the sum of 
two parts Ji and J2. Ji is proportional to the image part 
of the retarded Green function of the tip (or the density 
of state of the tip p ct (ui) = —^ImG ct ^(uj) ) and repre- 
sents the tunneling from the surface and absorbed atom 
into the tip. J 2 is proportional to the density of elec- 
trons (G^ t oc n) at the tip which represents the tunneling 
from the tip to the surface and absorbed atom. Experi- 
mentally, the shape of Kondo resonance in dl/dV curve, 
which is the mixing of Kondo resonance and others zero- 
bias structures, is not exactly the shape of Kondo reso- 
nance in the spectral function or density of state just like 
we have already discussed in the section B. The kondo 
resonance generally is modulated by a dip structure near 
zero-bias. 

We omit the Kondo Peak term Eq. J7J in the free Green 
function of the magnetic atom Eq. (JSJ) and calculate the 
dl/dV curve again. We can find that from Fig. the 
dl/dV curve is a dip. There is also a dip for empty sur- 
face when we remove the magnetic atom away from the 
surface. The dip generally becomes flat as increasing 
temperature. The behavior is similar to another impor- 
tant zero-bias phenomenon besides the Kondo resonance, 
the so-called zero-bias anomaly. [Hl| If we subtract the dip 
from total dl/dV curve we get the Kondo resonance Peak 
(solid- line). 

The orthogonality catastrophe [2!| of surface electrons 
prevents some many-body electronic transmission pro- 
cess. From the dl/dV curve we can find a local mini- 
mum near zero-bias, this is so-called zero-bias anomaly. 
The cotunneling mechanism has also contributions to the 
zero-bias anomaly even for nonmagnetic atom or quan- 
tum dot. j33, In this work the anomaly dip doesn't 
come from the orthogonality catastrophe because there 
aren't directly electronic interactions for the surface and 
tip. However the orthogonality catastrophe will con- 
tribute part of height of tunneling barrier from the sur- 
face to tip. The co-tunneling has been included in our 
model and self-consistent calculations. From Fig. ^ and 
the Feynman diagram in Fig.[3]for the Dyson equations, 
we can see that physics of co-tunneling is mainly included 
in the £3 self-energies for the tip, surface and magnetic 
atom, such as an electron jumps to spin-up state under 
the Fermi energy after a spin-up electron having already 
jumped to the surface. In order to prominently illustrate 
the effects of co-tunneling, we drop the S3 self-energies 
from the Dyson equations and calculate the dl/dV curve 
again. We find from Fig. [7| that the differential con- 
ductance decreases slightly. Remain dip structure is in- 
duced by subsequent single tunneling and higher power 
co-tunneling process embodied in the iterated calcula- 
tions of the Dyson equations. The co-tunneling events do 
happen with short time interval and they only play the 
role of virtually intermediate states in perturbed calcula- 
tions, but still have important influence on real physical 
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FIG. 8: The dl/dV curves change with (a) the amplitudes V 
from magnetic atom to tip, (b) B from surface to tip and (c) 
on-site Coulomb energy U at 40 K 



tunneling. 

In the absence of the Kondo resonance, no dip struc- 
ture near Fermi energies in the densities of state of the 
tip, magnetic atom and the surface show that the dip 
in dl/dV curve is dynamically tunneling phenomenon. 
Our results also indicate that the contributions of co- 
tunneling to zero-bias anomaly are not too large, but 
others tunneling events such as the subsequent single 
tunneling and the higher order co-tunneling contribute 
most parts of zero-bias anomaly. When Kondo resonance 
happens in the present of magnetic impurity, there is a 
dip structure in the density of state of surface conduct- 
ing electron [Fig|3(b)]. This contributes an additional 
dip structure to dl/dV curve, which mixes with the dip 
structure coming from dynamically tunneling. The dip 
structures have the widths with the same energy scale 
as the system temperature. If the Kondo temperature is 
higher than system temperature, the Kondo resonance is 
still modulated by the remain weaker dip structure [see 
Fig|7j and the Kondo resonance at first sight splits two 
peaks. The total shape of dl/dV curve is the mixing of 
two splitting peaks and main strong dip structure. 



E. Influence of tunneling amplitudes and electronic 
correlations 



The change of the tunneling current of STM is very 
significant when we change the distance from the tip to 
surface. The tunneling matrix V or B represents the tun- 
neling from the absorbed atom or the surface to the tip. 
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V and B are proportional to e~ d ^/ x , where d(r) is the 
distance from the tip to point r on the surface. The 
smaller is d the larger is V or B. Thus, when the tip con- 
tinuously accesses the surface or adsorbed atom, the pa- 
rameter V or B changes from smaller value to larger one. 
Fig. [Eta) illustrates the resonance peaks in dl/dV curves 
become stronger and asymmetric when V increases from 
small to large values. In calculations, we shift the Kondo 
resonance about SE — 0.005 eV above Fermi energy. We 
can also see that the Kondo resonance appears more 
prominent at strong tunneling in dl/dV curves. From 
Fig. |SJb) we also see that strong tunneling from surface 
to tip makes the shape of dl/dV curves are more sym- 
metric. This is because, at strong tunneling from the 
surface to tip, the role of the absorbed magnetic atom 
becomes less important. The electron-correlation inter- 
action plays an important role for Kondo effect. Based on 
the Yosida and Yamada's calculations EE3, for strong 
coulomb interaction (large U) the Kondo resonance sharp 
but weak with small area closed by the energy-axis. After 
having calculated the dl/dV curves, we find the Kondo 
resonance more peak alike but weak as having increased 
U from 2.0 eV to 4.0 eV. 



V. SUMMARY AND CONCLUSION 

The STM tunneling spectra are obtained using the 
non-equilibrium self-consistent calculations. The dl/dV 
curves have Fano's shapes and are close to the experi- 
mental tunneling spectra of 3d transition- metal atoms on 
novel-metal surface. The Kondo resonance in measure- 
ments dl/dV curves generally are the mixture of pure 
Kondo resonance and the dip structure due to others 
zero-bias anomaly such as the subsequent single tun- 
neling, cotunneling and other more complex tunneling 
events in the STM device. If Kondo resonance is strong 
the dl/dV curve is peak near zero-bias voltage and its 
shape is modulated by a valley or dip. Based on our 
model, the energy-window of tunneling electron has sig- 
nificant influence on the tunneling spectrum at very near 
zero-bias. The energy-window effect is unimportant for 
large bias voltage. 

On the one hand, Anderson model cann't capture all 
information in detail on the STM device, but it still cap- 
tures the most important properties of STM device. Our 
calculations only provide the main characters of Kondo 
resonance in STM device. More accurate calculations 
need to find a model more close to the real STM experi- 
ments and know the exact information on surface state. 

In this paper we have also developed a computational 
method of non-equilibrium electronic transport which is 
not only fit for STM device but also quantum dot, nano- 
molecule-device, hybridization system of normal metal 
and superconductor and others mesoscopic transport sys- 
tem. The self-consistent method can be extended to 
the system far beyond equilibrium by replace single- 
time integral and single- frequency Fourier transformation 



with double-times integral and double-frequency Fourier 
transformation. The method is suitable to the real time 
and real-space to study the pattern of tunneling spectrum 
on a surface. 
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APPENDIX A: NON-EQUILIBRIUM GREEN 
FUNCTION AND SCHWINGER-KELDYSH 
CONTOUR 

Schwinger and Keldysh [3^, Hj| introduced a time con- 
tour or closed time path shown in the schematic diagram 
below which includes a positive branch c+ from — oo to 
+oo and a negative branch c_from +oo to — oo. 

-oo c + 

' ) +oo 

-oo c - 

A new Green function is introduced on the contour, 
which is expressed as G c (t, t') = -i < T c A(t)^B(t') >, T c 
is an operator ordering the operators along the contour. 
If both t and t 1 are on the positive branch of the con- 
tour, T c is equivalent to the usual time-ordering operator 
and the Green function is the usual time Green function 
G T {t,t') = -i < TAffi B(t') >. If both t and t' are 
on the negative branch of the contour, T c is expressed as 
T which orders operators from +oo to — oo. The Green 
function expresses as G f (t,t') = -i < TA(t)^B(t') >. 
If t is on the positive branch and t' is on the negative 
branch, the order of operators is well organized and T c is 
omitted, the Green function G<(M') = i < A(t') ] B{t) > 
is the lesser Green function. If t is on the negative 
branch and t' is on the positive branch, the order of op- 
erator is converted and T c is omitted, the Green func- 
tion is the greater Green function and is written as 
G>(t,t') =-»< A(t)B{t'Y >. 

By introducing the time contour, Feynman rules, Wick 
theorem and the Dyson equation are the same as that at 
zero temperature. The only difference is that the time 
integration now is along the contour and the time can be 
at different branches of the contour. The Dyson equa- 
tion can be projected on the positive branch (— oo,+oo). 
If considering t and t' probably at different branches of 



10 



the contour, there are four Dyson equations which can regroup into a matrix equation |3J 
I 



QC 



G T -G< 



-9 



(Al) 



G c (t,t')=g c (t,t') + dh dt 2 g c {tM)^ c {tiM)G 2 {t 2 ,t')- 



(A2) 



Now the contour Green function can be calculated in 
the same manner as that at zero temperature. The only 
difference is that complex Green functions and self ener- 
gies are replaced by the complex matrixes of Green func- 
tions and self energies. In this work, we mainly study the 
electronic transport at small bias. We assume the STM 
work in the non-equilibrium stationary state. |3fij On the 
other hand the matrix Green functions define on only 
the positive branch. Thus, there is still the time trans- 
lation invariance after the system enters into the non- 
equilibrium stationary state. The matrix Green functions 
are dependent on the relative time interval G c (t — t'). In 
energy space, the Dyson equation is more simple under 
the non-equilibrium stationary state and written as 



G c (lu) = g c {uj) +g c {uj)Y 1 c {Lo)G c {uj). 



(A3) 



The retarded and advanced Green functions are calcu- 
lated using 



G» 



£ T H - E<(«) ; £ a H = £ r (w) - £> 
G t (lo) - G < (uj) = G>{lo) - G f(uj) 
G t (lo) - G>{oj) = G<H - G f (lj). 



In our calculations we have used following important 
properties of the matrix Green functions which can be 
derived from the symmetry of non-equilibrium Green 
function G<»(ui) = -G«»(w), G* T (u) = -G f {uj), 
G* t (lu) — G a (u>) and the Kubo-Martin-Schwinger condi- 
tion G<{uj) = -e-^ u) -^G > (cj) [3| 



= —2irif(uj)ImG r (uj), 

= -2i(n/(w) - l)/mG r (w) 

= ReG r {u) - i(2n/(w) - l)ImG r (uj), 

= -ReG r (uj) - i(2n/(w) - l)ImG r {tu) 



(A5) 



. where n. (cu) is the Fermi distribution function. 
APPENDIX B: DERIVATION OF THE 
TUNNELING CURRENT FORMULATION 

It is convenient to calculate the contour Green func- 
tion directly, not the lesser Green function. The mixing 
Green function in the current formulism (Eq[3J) can be de- 
fined on the contour and the full mixing contour Green 
function can be expressed using the Green functions of 
the tip, the adsorbed atom and the surface based on the 
methods in the reference. liOl 



(tl(t')d a ,(t)) c = \Y. J dt 0Gnn>s>At0,t')V:, ss ,g C «,(t ,t) 
n' ss' c 

(tl(t')p k<r ,(t)) c = ~ £ f dhG^Mta^Kk'^C'.J^ 



(Bl) 



s and s' are spin index. From the reference, ^3 if i n S inserted into the current formulation and the current 
A(t, t') — J dtiB(t, t%)C(ti, t') then formulation having been reorganized, the current is writ- 



A<{t,t')= / dt 1 B r {t 1 t 1 )C < {ti,t') + B<(t,t l )C a (t 1 ,t't$2) 

Jc+ 

After the contour Green function having projected 
onto the positive branch, the lesser Green function hav- 
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ten as (Jt) = J\ + J2, where 
J i = I^E f + dt G^ sa (t -t')T^<l(t ^m) 

nn' s<7 c 

= f^E Jmg^Imt^m] 



2e 



dwi^r[G d < r) (u)T ct < < >(a;)] 



2e 



2e 



du)ReTr[G cti<) (Lu)T ct< -^{uj)]. 



and the first equation converts to the second equation 
by using the Fourier transformations under the same 
assumption of non-equilibrium stationary state as the 
Dyson equations Eq. ( IA3(1 . The T matrixes are written 
as, 

hh' s' a' 

+ Ek^^C^^n 

s' a' 



q c p(<) anc | q cd (<) are th e f ree Green functions of the ab- 

^ kk' s' <r' ° V a' 

sorbed atom and the surface. In order to calculate con- 
tour Green function, at first, the self-energy is calculated. 
Based on Feynman diagrams [Fig[2], the self energies of 
the absorbed atom and the surface can be calculated us- 
ing the following formulations: 



n'n"a' a" 
■J'o'cj" 

y B , B* n „ g ct t (w) 

1" a' a" 

Ea^^^w+^m 

o-'tr" 

Z J maa> n aal a l ai V / 



k'k"a'o" 



n n a a 



a' 13' 



k'l'a'0' 



a and /3 are the spin index. The S3 self energies for the 
tip, adatom and surface are expressed as 



+ V A* V „, „B lRllR g ct „,» 5 C ?» 

= E 

k'l'ot'ff'a"0" 

+ Y B h , ,A* , „^ fl „„5 c f, , „M 

/ < mk'aa' V a " ™0 P 3 k ' I ' a ' ft 1 v ' ^ a" ff" y ' 

k'l'a'0'ct"/3" 

EaiflM = y A, ,.B* , „y*„ 3 c f, , ,H ff ct „ „(o;) 

mnfe'I'a'jfl'a"^" 

+ V V ,B „A* „ 3 ct , » 3 c f (u). 

/ ' mac' nk'p'a" l' fl" 13^ mna' P< v h' V a" $" V ' 

k'l'ot' j3'a" 0" 



(B6) 



which presents the cotunneling events in the STM device. 
A typical cotunneling event is that an electron jumps to 
the surface and leaved a hole on the magnetic atom, and 
subsequently a tip electron fills the hole on the magnetic 
atom [Fig. H an d Fig. Ej. For higher-order cotunneling 



processes include more electrons and tunneling events 
which generally are included in the iterated calculations. 
By compared the matrix expression of self-energy with 
the matrix expression of T matrix in the current formu- 
lation, T is just the lowest order approximation of the 



12 



exact self-energy. In our current formulation T replaces 
by the self-energy £ in Eq. I|B3|) . 

Ji = — / dujReTr[G ct{r HLo)^ ct(<) (Lu)} (B7) 
_o„ r+oo 

J 2 = — dLoReTr[G ct{<) (Lo)^ a \Lo)]. 

" i-oo 

When near zero bias the electronic tunneling is in classic 
forbidden regime, we use following current formula with 
energy cutoffs 

(Jt) = J1 + J2 (B8) 
Ji = ^ P + * duReTr[G ct( - r Hu)Z ct ^(uj)} 

J 2 = r +eVb ^i?eTr[G ct «)(w)S ct ( a )(w)]. 



where if positive bias voltage applies to surface V& > 0. 
/i is equal to /j,g the chemical potential of the surface. 
Because of numerical reason /zs sets to zero value to keep 
the position of Kondo peak unchange. The bias-voltage 
is applied by changing the chemical potential of the tip in 
our numerical calculations, [is keeps zero when applying 
a small bias voltage on the surface mean that applying 
a bias voltage to the tip V t =-Vb. 

In real calculations, we add a minus ahead the current 
formulas Eq.®, Ea. (|B7l) and Ea. (|B8l) ) to make the cur- 
rent Jt = J\ + J2 > at positive bias once the lower 
bound and up-bound of the integrates are clearly writ- 
ten in the current formulas. We can adopt the chemical 
potential langue just like in Ref |39( without introduc- 
ing additional minus in the current formulas. The basic 
reason of the difference is that the increase of chemical 
potential 5fi of electron is equal to -eSVb 
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